Oscillating Quintom and the Recurrent Universe 
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Current observations seem to mildly favor an evolving dark energy with the equation of state 
getting across -1. This form of dark energy, dubbed Quintom, is studied phenomenologically in this 
paper with an oscillating equation of state. We find oscillating Quintom can unify the early inflation 
and current acceleration of the universe, leading to the oscillations of the Hubble constant and a 
recurring universe. Our oscillating Quintom would not lead to a big crunch nor big rip. The scale 
factor keeps increasing from one period to another and leads naturally to a highly flat universe. The 
universe in this model recurs itself and we are only staying among one of the epochs, in which sense 
the coincidence problem is reconciled. 

PACS numbers: 98.80.Cq 



In 1998, two groups Q, Q independently showed the 
accelerating expansion of our universe, which established 
the existence of dark energy where the equation of state 
is less than —1/3 today. The simplest form of dark en- 
ergy is the cosmological constant. A cosmological con- 
stant which leads to current acceleration would encounter 
many theoretical problems, such as the fine-tuning prob- 
lem and the coincidence problem [3j. A light scalar field 
of quintessence 0, U @ which evolves with time is to 
some extent likely to resolve the coincidence problem. 
The model of phantom|7j has also been put forward 
which leads to an equation of state w < — 1. Normally 
quintessence with a canonic kinetic term can only have 
w > — 1, meanwhile k-essence.S) can have both w > — 1 
and w < — 1. For the cosmological constant and many 
quintessence models the event horizon would lead to a 
potential incompatibility with the string theory, mean- 
while for models where w < — 1 one would get the Big 
Rip 9| of the universe. 

The accumulation of the current observational data has 
opened a robust window for probing the recent behavior 
of dark energy. Measurements from type la Supernova 
(SNe la), the Cosmic Microwave Background (CMB) Ra- 
diation, Large Scale Structure (LSS), weak lensing, clus- 
ters and galaxies all contain the imprints of dark energy. 
Specifically the recently released first year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) measurement [T(i| . 
the Sloan Digital Sky Survey (SDSS) measurement of 
the three-dimensional power spectrum|ll| and most im- 
portantly, the recent discovery of 16 SNe la E2 with the 
Hubble Space Telescope during the GOODS ACS Trea- 
sury survey, together with former SNe la data have pro- 
vided the most precise up-to-date measurements of dark 
energy. Many other sources have also been studied to 
make crosschecks and reveal new physics on the dark en- 
ergy and the concordance cosmological model. Recently 
many authors in the literature have fitted the behavior 
of dark energy using various parameterizations. The best 
fit value of the equation of state is found to be less than 
— 1 when using SNe la 

[3 El 13 El 13 113 or the X-ray 
mass fraction data^^ with a cosmological constant still 



well within the central regions 1 . It has also been pointed 
out by many authors that an evolving dark energy is in- 
deed more favored than that with a constant equation of 
state. With the simplest linear parametrization of the 
equation of state w: 



w{z) — wq + w'z, 



(1) 



SNe la alone favors a w larger than — 1 in the recent past 
and less than —1 today, disregard using the prior of a 
fiat universe [3 El El El El or not [l7|. When other 
"longer-armed" observations, e.g. CMB and LSS data, 
are also taken into account, one often takes the risk of 
poor parametrizations. Taking the linear parametriza- 
tion as a concrete example, it gives a good approxi- 
mation for the late behavior of the dark energy. How- 
ever, at large redshift the amplitude of w increases too 
much and hence CMB and LSS data would restrict w' 
to be near zero0,|2(j> spoiling the recent behavior of 
dark energy. Such a correlation also happens for other 
parametrizations. The most proper way of measuring 
the recent behavior of the dark energy is the principal 
component analysis |Elj. Huterer and CoorayEjj based 
on this method and developed a simple uncorrelated and 
nearly model-independent band power estimates of w and 
its density as a function of the redshift. By fitting to the 
recent SNe la data they found a marginal (2-cr) evidence 
for w(z) < — 1 at z < 0.2. Even though SNe la data are 
still currently limited, such a method is potentially very 
powerful. The parametrization of the dark energy re- 
mains a paradox today, especially when confronted with 
observations |22| . However, in any case if a time- varying 
dark energy model significantly reduces the % 2 value and 
is statistically significant, it gives strong implications for 
dark energy metamorphosis. In Refs. [23ll24| the authors 
use their physical parametrizations of the dark energy 
and fit the models to the fully CMB, LSS and SNe la data 



1 We get the similar results when using the recent data of radio 
galaxies 19] as well as the gold samples of SNe la data ll2l 
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set, they also find that an evolving dark energy is favored 
indeed with w(z) < —1 today and w(z) > —1 in the re- 
cent past. However, one would have to introduce several 
additional parameters to give a physical interpretation 
for the behavior of the dark energy, resulting into a very 
low statistical significance of an evolving dark energy. 
Most of the works in the literature have not taken into 
account the possibly non-negligible contribution of weak 
lensjng[25| to SNe la. When the flux averaging method 
|26l |27| has been applied the preference of a nonzero w' 
is much reduced. In Ref.^U we have also applied the 
linear parametrization to the X-ray mass fraction data 
only and find that a negative w' is favored. 

If the requirement of an evolving dark energy with w 
getting across — 1 still holds on with the accumulation of 
observational data, this would be a great challenge to the 
current cosmology. In the quintessence model, the equa- 
tion of state w = p/ p is always in the range — 1 < w < 1 
for V(Q) > 0. Meanwhile for the phantom which has 
the opposite sign of the kinetic term compared with the 
quintessence in the Lagrangian, one always has w < — 1. 
Neither the quintessence nor the phantom alone can ful- 
fill the transition from w > — 1 to w < - 1 and vice versa. 
Although for k-essence@ one can have both w > — 1 and 
w < —1, it has been lately considered by Ref|2^ that it 
is very difficult for k-essence to get w across —1 during 
evolving. In Ref. |l5j we proposed a new scenario of dark 
energy dubbed Quintom. Quintom describes those forms 
of dark energy which can get w cross — 1 during the time 
evolution. 

When dark energy is not merely the cosmological con- 
stant it gives rise to many interesting consequences. 
There are possibly some connections between the early 
inflation and the current acceleration of our universe. 
Peebles and Vilenkin have put forward a model of 
quintessential inflation to unify the two epochs [2j| . Do- 
delson, Kaplinghat and Stewart gave a simple form of 
oscillating quintessence which can provide a natural so- 
lution to the coincidence problem 30] . In this paper we 
propose a phenomenological model of oscillating Quin- 
tom and study its implications in cosmology. We will 
show that in our model two accelerating epochs are re- 
lated and the coincidence problem can be alleviated. 

The phenomenological Quintom model we consider in 
this paper takes the following equation of state, 

wQxia) = wo + w\ cos[Aln(a/a c )], (2) 

where a is the scale factor, u>o, Wi, A and a c are param- 
eters which we will take specific values for the detailed 
study below. If for a sufficiently large period the equation 
of state remains unchanged during the late time evolving 
it would fit all current data well for w < — 1. Defining 
X(lna) = px{lna)/px(0) where the subscript X stands 
for Quintom and '0' for today, one gets: 

X(\n a) = a - 3{1+Wo) exp{^i [sia(A In — ) + sin{A In a c 

A a c 

(3) 



For a numerical discussion we fix the model parameters 
and consider a specific w given below 

47T 

w(hxa) = -1- 1.5cos(0.0321na ). (4) 
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In this model the w at present time is w = —1.26 ( we 
have checked that if tuning a c to get w more closer to 
— 1 the picture on the evolution of the universe will not 
be changed from our current choice of the parameters), 
which is well within the limit of current observations. In 
Fig.l we show the evolution of the universe filled with 
only the matter of the oscillating Quintom. As we've 
fixed wq = — 1 , one can see from Eq|3| that the energy 
density, hence the Hubble parameter H simply oscillates 
with a long period. Since the time is a simple integration 
of dhia/H and X(lna) evolves periodically with lna, 
time increases equally for every span of lna. One can 
also easily find from Fig.l that as X is an integration of 
w(lna) they differ in the phase of oscillation. Since the 
Hubble parameter increases gradually during phantom- 
dominated phase of dark energy, this will make the uni- 
verse evolve into a regime with high energy again 2 , which 
may be regarded as the recurrence of an early phantom 
inflation [3l| of our universe, (see [12 for a earlier study) 
where the initial perturbation in the horizon will exit the 
horizon, and reenter the horizon after the transition to a 
late-time expanding phase with the ordinary matter con- 
tents. Thus our model unifies the inflation and dark en- 
ergy into a single process, in which the universe evolves 
periodically. While we live in one period, the present 
acceleration with w < —1 is just a start of the phan- 
tom inflation for the next period of universe. During this 
process, all ordinary forms of matter and radiation are di- 
luted by inflation. But in the next period new structures 
will be formed again by the fluctuations generated during 
the pre sent acceleration, i.e. phantom inflation [3lJ (see 
|33l l34l ] for a detailed investigation of primordial pertur- 
bations from phantom phase). The universe can also be 
reheated with help of the reheating mechanisms such as 
instant preheating)^ and curvaton reheating|36|. 

We show now that our model works the same as de- 
scribed above when including matter and radiation. As- 
suming f2 TO = 0.3 and fix — 0.7 today with the Hubble 
parameter h = 0.7 ( this determines the abundance of 
radiation ) and a negligible curvature, we show in Fig. 2 
the evolution of the universe for 1-period. Here we as- 
sume that the radiation starts to dominate our universe 
promptly at the moment when ln(a/ao) = —60 ( ao is 
set to be unity today ). From Fig. 2 one can see that 
firstly the universe is just like today with fix = 0.7 
and w = —1.26, with time w evolves downward, px in- 
creases and p m decreases, then the universe begins in- 
flating. When lna ~ —150, w starts to increase. And 



_ 2 This owes to the special features of phantom, in which the null 
energy condition is violated. For the oscillation of state equation 
without vi < — 1, this case does not occur. 
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FIG. 1: Cosmological consequences with only an oscillating 
Quintom it; (In a) = -1 - 1.5cos(0.0321na - 



at In a 



-90 we have w 



-i, inflation stops and the 
energy density continues dropping. At In a = —60 the 
universe enters radiation domination epoch, px starts 
to increase recently and begins to dominate the universe 
with fix = 0.7 today. The universe starts the inflation 
since today and is recurring into the next period. 



multi periods. For each period we assume that the ra- 
diation starts domination around the same time. As ra- 
diation and Quintom domination epoch contribute little 
to time compared with when matter is non-negligible, 
one can see from Fig. 3 that time only increases equally 
around every short epoch of matter domination. 
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FIG. 3: Cosmological consequences with an oscillating Quin- 
tom io(lna) = -1 - 1.5cos(0.0321na - ^f) and fi m = 0.3, 
tlx = 0.7 and h = 0.7 today, within several periods. 
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FIG. 2: Cosmological consequences with an oscillating Quin- 
tom iu(lna) = -1 - 1.5cos(0.0321na - ^) and tt m = 0.3, 
Qx = 0.7 and h = 0.7 today, within only one period. 

In Fig. 3 we show the evolution of the universe with 



We should point out that for our phenomenological 
parametrization in Eq. (J2J) and Eq.Q) to realize our re- 
current picture of the universe, the value of wq has to 
be extremely close to -1, otherwise the next and previous 
"cycles" would be different from the current one and we 
would also be encountered with too high or low inflation 
energy scales, as can be clearly seen in Eq.©- While for 
the other several parameters less fine tunings are needed: 
the period is determined by A and the variation of the 
energy density is by loi and A, and finally, the value of 
a c somewhat determines the present value of w. To real- 
ize the recurrence picture one requires that inflationary 
epoch is satisfactory and hence we have 2ir/A ~ 100, on 
the other hand we need Hi n f/Mpi < 10~ 5 to satisfy the 
primordial gravitational wave constraints, these will put 
considerable constraints on w\ . For our dynamical model 
of dark energy as the period of oscillation is required to 
be quite large, the crossing of the cosmological bound- 
ary appears at some high redshifts where dark energy is 
typically negligible compared with the background mat- 
ter components, where such a crossing behavior is hard 
to be detected by the observations. Nevertheless such 
a crossing is crucial to realize our picture of oscillating 
quintom. The dynamics of dark energy is not well im- 
printed on the redshift range covered by current SNe la 
surveys, where a running of the equation of state dw/dz 
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would be typically of order 0.01, this is different from 
the model of Quintom investigated earlier where a signif- 
icant dynamical dark energy behavior displays in a very 
low and short redshift range [l5|. For the specified exam- 
ple in Eq.||IJ the running of w is ~ 0.02, which is difficult 
to be excluded by future observations like SNAPI Sflillfl . 
Joint Efficient Dark-energy Investigation ( JEDI) |4lL l42j 
or Dark Energy Survey (DES)|43j. On the other hand, 
if a significant dynamical behavior of dark energy were 
verified by future observations, this would be a strong 
argument to rule out our model. Again we should point 
out that Eq. is one of the typical examples where the 
current equation of state is quite adjustable which does 
not affect our picture. To satisfy the current observations 
we need the equation of state close to —1 today, which 
in turn requires in Eq.Q the term w\ cos [A In a c ] to be 
close to zero. 

During the phantom-dominated phase, the energy den- 
sity increases as well as the Hubble parameter. If there is 
not a transition, the exit to observable universe is hardly 
possible, which is also a problem puzzling phantom infla- 
tion 3 . But in our model it is avoided by the oscillation 
of w. In some se nse, our model is similar to the Pre 
Big Bang scenario |44| ( see 145] for a review) and the re- 
cently proposed cyclic model [4|| where the primordial 
perturbations are seeded by the contraction phase with 

3 The exit may be implemented by using growing wormholes and 
ringholes l37l . which can be regarded as a connection bet ween 
the phantom phase and late-time expanding phase, see l38l for a 



the increase of the energy density. To exit and enter late- 
time observable universe, a non singular bounce has to 
be required. But different from those which lead to the 
big crunch of the universe, in our model the universe is 
always in the expanding phase with periodicity. In this 
case the "bounce" actually means a transition from one 
expanding phase to another. 

In summary we have proposed in this paper a phe- 
nomenological model of Quintom with an oscillating 
equation of state. Such a model naturally unifies the 
early inflation and the current acceleration of the uni- 
verse, leading to oscillations of the Hubble constant and 
a recurring universe. The model of oscillating Quintom 
does not lead to a big crunch nor big rip. The universe 
just recurs itself with the scale factor increasing always 
and we are only staying among one of the epochs, in 
which sense the coincidence problem is somewhat recon- 
ciled and fife is predicted to be highly close to zero 01 • 
Further studies on model building, perturbation and re- 
heating are currently under investigation^^. 
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